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Abstract 

We describe a new experimental setup for the detection of magnetic circular dichro- 
ism with fast electrons (emcd). As compared to earlier findings the signal is an order 
of magnitude higher, while the probed area could be significantly reduced, allowing 
a spatial resolution of the order of 30 nm. A simplified analysis of the experimental 
results is based on the decomposition of the Mixed Dynamic Form Factor S(q, <f , E) 
into a real part related to the scalar product and an imaginary part related to the 
vector product of the scattering vectors q and if. Following the recent detection of 
chiral electronic transitions in the electron microscope the present experiment is a 
crucial demonstration of the potential of EMCD for nanoscale investigations. 



1 Introduction 

The observation of circular dichroism with electron probes has been considered 
impossible except with spin polarized electron probes. In 2003 it was suggested 
that this was not the case [lj. The magnetic transitions that give rise to X-ray 
magnetic circular dichroism (xmcd) contribute to the imaginary part of the 
Mixed Dynamic Form Factor (mdff) [2j for inelastic electron scattering. Since 
this quantity can be measured in the transmission electron microscope (tem) 
under particular scattering conditions, we predicted that detection of XMCD 
should be feasible in the TEM. We called the predicted effect Energy-Loss 
Magnetic Chiral Dichroism (emcd). 
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In the first conclusive experimental demonstration of EMCD on Fe [3 J it was dis- 
covered that the effect is smaller than XMCD. In that experiment the dichroic 
signal was close to the noise threshold in the then chosen geometry and the 
area of investigation was approximately 100 nm in radius. 

Here we present a new experimental setup that enhances the count rate by 
an order of magnitude and reduces the probed area by another factor of five, 
thus opening the way to applications of EMCD on the nanometric scale. 

The MDFF S(q, <f , E) is the essential quantity describing EMCD. It has been 
used in the description of interference of inelastically scattered electrons (e. g. [2ll4ll5] ). 



2 The Mixed Dynamic Form Factor (MDFF) 



The semi-relativistic double differential scattering cross section for inelastic 
electron scattering (DDSCS) in the plane wave Born approximation is [6j 
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where a is the Bohr radius, ki (kf) is the wave number of the incident (outgo- 
ing) probe electron, q = — kf is the wave vector transfer in the interaction 
and E the energy loss. S(q, E) is the dynamic form factor (dff). 

Interference between inelastically scattered electrons in the diffraction pattern 
will occur when the probe electron consists of two or more mutually coherent 
plane waves [HIT] . Experimentally, this can be realized by a biprism [8] or by 
any other beam splitter. It was shown experimentally that the crystal itself can 
be used as a beam splitter for inelastic electron scattering [9j. In the crystal the 
probe electron is a superposition of Bloch waves which, in turn, are coherent 
superpositions of plane waves defined by the allowed Bragg reflections. 

For the sake of clarity, we consider the simplest case here, namely the su- 
perposition of two plane waves with complex amplitudes respectively. 
Technically, this situation is approximated in electron diffraction by the two- 
beam case, the most important plane waves being the incident one and a single 
Bragg scattered wave. 

The DDSCS is then [2] 
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Here, q = ki — kf, q' = k[ — kf are the wave vector transfers from the two 
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incident plane waves k i) k[ to kf. Since i) the 2p — > 3d transitions are the dom- 
inant ones due to the shape of the density of states and ii) due to the localized 
character of < z|, the matrix elements (i\e l ^ R \f) are mostly determined by an 
area within small R values (compared to lattice parameters), we can use the 
dipole approximation for the MDFF [4] 

S(q, q', E) dip = ^\^\f)(f\^R\i)S{E + E t - E f ). (3) 

if 

R is the 3-vector operator (ri,r2,r3) of the one-electron scatterer with initial 
and final wave functions |/). 

Eq. [2] consists of two direct terms, each resembling the angular scattering 
distributions centered at the incident and the Bragg scattered plane wave 
directions, and an interference term. Eq. [2] is formally equivalent to the ex- 
pression for intensity in the double slit experiment. It should be noted that 
the diagonal element of the MDFF is the DFF, S(q, q. E) = S(q, E). 

The MDFF describes the mutual coherence of transitions with energy transfer 
E and momentum transfer hq, Tiq 1 [4J (Two different momentum transfers can 
occur in one transition with finite probability when the incident or the outgo- 
ing electron is not a single plane wave. In the present case Ai\ki) + A 2 \k f j ) is 
such a basis function, and the measurement collapses the probe electron into 

With the matrix elements 

rjk = E<*N/></N^(£ + E f ). (4) 

if 

of the transition matrix R = {rjk} the MDFF, eq. [i]can be written as 

S{q,q',E) dip = qRq' . (5) 

For isotropic systems the transition matrix degenerates to a quantity pro- 
portional to the unity matrix [2j. This case was discussed in the context of 
ionisation fine structure and dynamical diffraction [TOUT] . 

Anisotropy can be induced by a lattice of lower than cubic symmetry, or by 
magnetic fields. These can be internal or external, then speaking of natural or 
magnetic dichroism jllj . It is well known that with photon scattering linear as 

1 It should be noted that collapsing the probe electron into \kf) does not exclude 
the possibility of interference between outgoing beams with different wave vectors; 
any outgoing beam that is Bragg scattered to \kf) before leaving the crystal can 
produce interference detectable by the setup described here. 
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well as circular dichroism can be measured. This technique is largely applied 
with external magnetic fields. Linear magnetic dichroism shows up as an uni- 
axial anisotropy and can be measured with angle resolved inelastic electron 
scattering, tuning q parallel or perpendicular to the anisotropy axis. This is 
equivalent |12] to the tuning of linear polarization of the photon in XANES 
experiments. 

It has been thought that circular magnetic dichroism cannot be detected with 
electrons except with spin polarized ones. But let us recall that in XANES the 
photon does not couple directly to the spin of electrons but to the angular 
momentum of the excited atom, and the effect becomes visible by the spin- 
orbit coupling [llj. So there is no reason that spin polarized electrons are 
needed for detection of circular dichroism in electron energy loss spectrometry 
(eels). Rather, in the inelastic electron interaction that is equivalent to an 
XMCD experiment, the virtual photon that is exchanged must be circularly 
polarized. 

The MDFF, eq. [5] can be written in a different form when we specify the 
magnetic field direction as the positive r 3 axis p~| and write q = (q±,qs). A 
little algebra shows that 

1 i 

S(q,q , ,E) dip = -(r + ++r__)q ± 'q ± , + r 00 q 3 q , 3 + -(Y+ + -r__) (q ± x q± f )-e 3 (6) 

where e 3 is the unit vector in direction of the r 3 axis, and we have used the 
transition matrix elements in terms of the spherical components i? +5 _ ? o of the 
3-space operator 

r ++ = J2(i\R+\f)(f\R+\i)5(E + E f ). (7) 

if 

and similar for all other combinations. They relate to the Cartesian com- 
ponents by the transformation rules for vector spherical harmonics [13j. All 
matrix elements and vector components are real in eq. [6j In this form we have 
separated the MDFF into a real component proportional to the scalar prod- 
uct of the wave vector transfers gl, q'j_ and an imaginary part proportional 
to their vector product. This structure is equivalent to that of the polariza- 
tion tensor used in XMCD, which decomposes into a scalar part (uneffective 
in dichroic experiments), a second-rank irreducible part detectable by linear 
dichroism and a pseudovector part sensitive to magnetic moments [TT1fT4] and 
reference therein. 

This form of the MDFF allows a description of the scattering geometry for 
EMCD detection in the TEM. In passing we note that the imaginary part van- 
ishes if the magnetic transitions (Am = ±1) are degenerate. Only when the 

2 In the TEM this is usually also the optical axis. 



4 



presence of a magnetic field in r% direction lifts the m-degeneracy will we see 
an effect. For a transition with fixed energy loss E the operators R+ and R_ 
will then contribute with different oscillator strengths, and the MDFF eq.[6]will 
acquire an imaginary part. Its sign depends on which transitions are allowed 
by the selection rules. 

The imaginary part can be interpreted as the difference in probability to 
change the magnetic quantum number by ±1. It thus describes the differ- 
ence in response of the system to left- respectively right-handed circularly 
polarized electromagnetic fields p~[ 

The scattering vector q± in the diffraction plane is perpendicular to the mag- 
netic field vector. We assumed already that the magnetic moments of the 
scatterer are aligned parallel to the optical axis r 3 in the strong magnetic field 
2 T) of the objective lens of the microscope. We can now evaluate the 
DDSCS for specimens showing magnetic circular dichroism. 

If in eq. [2] we write the phase shift (j) explicitly as A1A2 = | -Ai 1 1 -A2 1 • e - ^, 
inspection of eq. [6] then reveals that a phase shift (j) ^ nir between the two 
incident plane waves is needed in order to activate the imaginary part of the 
MDFF. A phase shift of ±7r/2 is recommended since in this case the real part of 
the MDFF disappears in eq. [2j and only the imaginary part survives. In a two- 
beam case with such a phase shift the pseudovector part contributes with its 
full magnitude and gives rise to an asymmetry in the scattering cross section 
of a magnetic transition such as the L edges of the ferromagnetic d-metals. 



3 Experiments and simulations 



We performed our experiment on a Co single crystal electropolished sample 
with a FEI Tecnai F20-FEGTEM S-Twin operating at 200 keV and equipped 
with a Gatan imaging filter (GIF). The dichroic signal is obtained by first 
tilting out of the [0001] zone axis to a two-beam case where only the 0000 and 
1010 reflections are strongly excited. Following the procedure illustrated in [3], 
a selected area aperture (SAA) is used to delimit a region of about 100 nm 
radius and 18 ± 3 nm thickness. The corresponding diffraction pattern is then 
projected onto the 2 mm spectrometer entrance aperture (SEA). Drawing a 
circle with a diameter of G with the 0- and G-reflections on the left and right 
side respectively, the strongest dichroic signal can be expected at the top and 



3 An imaginary part of the MDFF signifies that time inversion symmetry is broken. 
In fact this symmetry breaking relates to the angular momentum operator. Under 
time inversion its direction is reversed. In the presence of a magnetic field this is no 
longer a symmetry operation. 
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Fig. 1. Left: the electron beam creates a diffraction pattern in the focal plane of 
the objective lens. Placing the SEA in the diffraction plane will select one direction 
for the outgoing electron beam. The red circles indicate positions A and B of our 
experimental setup. Right: in reciprocal space the Ewald sphere is defined by the 
incident wave vector lq. If a reciprocal lattice vector (g) lies on the surface of the 
sphere, the Bragg spot will be strongly excited (as is the case for kj in the figure). 
The projection of the sphere on the diffraction plane is called Laue Circle Centre 
(LCC, equal to g/2 in the figure). 

bottom points A and B of this Thales circle (fig. [T]) where the scalar product 
q • q ' is zero and the pseudovector qxq f maximises the imaginary part of the 
MDFF. The camera lenght and projection coils are then adjusted so that the 
SEA will be centered at these two points (first A, then B) in the Thales circle 
and determine a collection angle of 2-4 mrad. Two spectra are then acquired 
sequentially (fig. [2]). With an acquisition time of 60 sec per spectrum and an 
energy dispersion of 0.5 eV/channel the intensity at the L 3 peak was about 
2,500 counts (after background removal). 

In order to improve the signal-to-noise ratio a new experimental setup was 
devised. As detailed above, the sample is tilted out of the [0001] zone axis 
to a two-beam case where only the 0000 and 1010 reflections are observed in 
the the energy filtered diffraction pattern, which is then projected onto the 
spectrometer entrance aperture (SEA). Using a rotational sample holder, the 
reciprocal lattice vector G is then aligned parallel to the energy dispersive 
axis of the CCD camera, so that a q-E diagram can be recorded as depicted in 
fig. [3| The quadrupoles of the energy filter collapse (integrating the signal in 
the q x dimension) the circular area to a line in q y when the system is switched 
to spectroscopy mode. This first modification allows us to record not only 
both spectra A and B with a single acquisition, but the entire range of spectra 
with different q y values comprised within the 2 mm SEA. It should be noted 
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Fig. 2. EELS spectra of the Co £2,3 edge showing EMCD using the detector shift 
method. 

however that the integration area in the q x dimension is different for every 
value of \q y \. 

A second modification consists in a different method [15] to obtain a spot 
like inelastic diffraction pattern: the beam (with a convergence semi-angle of 
a = 2 mrad) is focussed onto a 18 ± 3 nm thick area of the Co specimen 
which is then shifted upwards from the eucentric position by z = 9.25 /im. 
The diameter of the illuminated area is d = 2az = 37 nm, accurate to 5%, in 
the present experiment (fig. [i]). 

Preliminary experiments show that with smaller z-shifts the illuminated area 
can be reduced to less than 10 nm radius, at the moment with untenable 
distortions of the diffraction pattern. With a C s corrector or a monochromator 
a spatial resolution of 10 nm or less should be attainable. 

Similarly to XMCD we define Aa as the difference between spectra with op- 
posite helicity and a as their average. The dichroic signal is then the relative 
difference Aa/a in the scattering cross section when the sign of the pseu- 
dovector part changes Q This is obtained by tracing the spectral intensity at 
points A and B in fig. |3| and taking their difference, as shown in fig. [5| With 
an acquisition time of 15 sec and an energy dispersion of 0.3 eV/channel the 
intensity at the L 3 peak was about 13,500 counts (after background removal). 

The modified scattering geometry provides a count rate per eV which is an 
order of magnitude higher than the one achieved in the previous configura- 
tion [3j, thus improving significantly the signal to noise ratio. This is essentially 
caused by the fact that when no SAA is used and the beam is focussed only on 
the area of interest, all the electrons emitted from the gun contribute to the 

4 Another common definition of the dichroic signal is the ratio between the difference 
and the sum of spectra with different helicity, which is a factor of 2 smaller than the 
one used here. 



7 




energy loss [eV] 



Fig. 3. Left: Energy filtered diffraction pattern at 784 eV energy loss using a slit 
width of 20 eV. The sample was oriented in the two-beam case, capturing both 
reflections within the 2 mm SEA. Right: q-E diagram of the Co edge in chiral 
conditions (for details see text). 



Specimen 1 


Jr— Illuminated Area 










z-shift 



Image Plane 



Fig. 4. Geometry of the experiment: When the specimen is shifted upward from the 
eucentric position (Image Plane) by z the illuminated area is a disk with radius az. 
Bragg scattering will cause a diffraction pattern to appear in the Image Plane. 

detected signal. In the formerly used geometry a nearly parallel incident bun- 
dle illuminated a large area of the sample of which only a small fraction could 
be used. This effectively reduced the intensity by which the area of interest is 
illuminated, i.e. a large part of the incident electrons did not contribute to the 
signal. The increase in the count rate per eV allows us to reduce the acquisi- 
tion time, thus limiting the effects of beam instability, specimen and energy 
drift. The shorter acquisition time, combined with the finer energy dispersion, 
improves the energy resolution with which the edges are recorded. In the 
older setup the L 3 has a FWHM of 7 eV (fig. [2]), compared to the 3.6 eV 
achieved with the new method (fig. [5]). 
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Fig. 5. EELS spectra extracted from the q-E diagram (fig. [3]) by line traces parallel 
to the energy axis and passing by the positions A and B. The difference is the 
dichroic signal stemming from ^s(S(q, q\ E)). 

Ab initio DFT simulations of the dichroic signal were performed using an 
extension [16] of the WIEN2k package [17] developped for this purpose. In the 
simulation the effects of thickness, tilt of the incident beam, position of the 
detector were included as well as the integration over q x in the range dictated 
by the use of a circular SEA. Up to 8 beams were used for the calculations 
of the MDFFs. A comparison with the experiment is given in fig. [6] for the 
L 3 edge of Cobalt. The agreement is very good between -0.8 and 0.8 G with 
some discrepancy appearing at larger scattering angles. This can be due to 
the faint Bragg spots outside the systematic row (which are neglected in the 
simulations) and to the fact that the SEA is not exactly in the spectral plane 
of the energy filter. The error bars correspond to the (2a) Poissonian noise 
calculated for the theoretical signal using the number of electrons contributing 
to the signal as determined from the experimental data. 



4 Conclusions 

The strong chiral effect observed in the Co L 2 ,3 edge shows that EMCD can be 
measured with reasonable collection time in the TEM. The convergence angle 
of 2 mrad is obviously not detrimental for the necessary constant phase shift 
between the unscattered and the Bragg scattered electron waves. The z-shift 
of the specimen allows to control the illuminated area, and with optimized 
conditions a lower limit of < 10 nm appears realistic. This would define the 
lateral resolution in scanning mode. In order to achieve this goal some techni- 
cal problems such as the stability of the beam in the magnetic field, constant 
z-shift or decoupling of the scan coordinate from the positioning of the diffrac- 
tion pattern must be solved. The simple concept described above should be an 
incentive for novel dichroic experiments in the TEM. It also shows that EMCD 
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POS B 
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Fig. 6. Dichroic signal at the Cobalt L3 edge as function of the scattering angle 
q (in unit of G) in the direction perpendicular to the Bragg scattering vector G. 
Positions A and B (see fig. [3]) correspond to points q/G = 0.5 and -0,5 respectively. 
Comparison between experimental data (triangles) and ab initio simulation. The 
error bars correspond to simulated Poissonian noise (2a). 

can be complementary and competitive with traditional or new [18j XMCD 
techniques. 
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